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Brief  title:  Correlation  coefficient  images 

Abstract: 

In  ultrasonics,  image  formation  and  detection  are  generally  based  on  signal  amplitude.  In  this 
paper,  we  describe  an  amplitude  independent  approach  for  image  formation  and  detection  based 
on  the  similarity  of  adjacent  signals.  Signal  similarity  is  quantified  in  terms  of  the  correlation 
coefficient  calculated  between  A-scans  digitized  at  adjacent  measurement  positions.  Correlation 
coefficient  images  are  introduced  for  visualizing  the  similarity  in  measured  A-scans.  In 
backscatter,  the  approach  reveals  defect  signals  buried  in  noise  by  showing  regions  of  increased 
correlation.  In  pitch-catch  or  thru-transmission,  the  approach  reveals  defects  by  showing  regions 
of  decreased  correlation  due  to  signal  distortion  caused  by  interaction  of  the  beam  field  with  the 
defect.  Correlation  coefficient  and  C-scan  images  are  shown  to  demonstrate  flat-bottom-hole 
detection  in  a  stainless  steel  annular  ring  and  crack  detection  in  an  aluminum  aircraft  structure. 
Simulated  data  are  used  to  show  the  detection  of  planar  defects  at  very  low  signal-to-noise  ratio. 

PACS  numbers:  43.20.Hq,  43.20.Ye,  43.20.-f,  43.35.Cg 
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I.  Introduction 


Typical  ultrasonic  detection  in  nondestructive  evaluation  (NDE)  is  amplitude  based.  In 
the  simplest  case,  an  echo  signal  above  a  threshold  or  a  transmitted  signal  below  a  threshold 
form  the  basis  for  detection.  C-scan  images  are  used  to  display  peak-amplitudes  throughout  a 
scan.  Based  on  a  C-scan  image,  regions  of  high  or  low  amplitude  can  be  identified  either 
manually  or  using  a  computer  decision  based  approach.  A  weakness  of  amplitude  based 
approaches  is  their  direct  and  inherent  sensitivity  to  non-defect  related  amplitude  changes 
associated  with  the  measurement  system  or  sample.  In  this  paper,  we  focus  on  a  complimentary 
detection  approach  which  relies  on  the  correlations  between  adjacent  A-scans  as  a  basis  for 
detection. [1-3]  Correlation  coefficient  images  are  introduced  for  visualizing  the  similarity  in 
measured  A-scans.  Rather  than  using  signal  amplitude,  this  approach  bases  detection  on  the 
similarity,  or  lack  of  similarity,  between  adjacent  signals.  In  backscatter,  this  approach  reveals 
defect  signals  buried  in  noise  by  showing  regions  of  increased  correlation.  In  pitch-catch  or  thru- 
transmission,  the  approach  reveals  defects  by  showing  regions  of  decreased  correlation  due  to 
signal  distortion  caused  by  interaction  of  the  beam  field  with  the  defect. 

The  basic  approach  for  forming  correlation  images  as  described  in  this  paper  uses  the 
same  A-scans  which  form  the  basis  of  C-scan  images.  In  either  case,  the  process  begins  with 
measured  A-scans  written  into  a  three-dimensional  matrix  with  the  rows  and  columns  of  the 
matrix  registered  with  the  measurement  locations  and  with  the  A-scans  written  into  the  third 
dimension  (see  Fig.  1).  A  C-scan  image  is  typically  formed  by  establishing  a  time  gate,  finding 
the  maximum  absolute  value  (or  maximum  min-to-max  deviation)  within  the  time  gate  for  each 
A-scan,  writing  these  maximum  values  into  a  two-dimensional  matrix  which  mimics  the  scan 
pattern,  and  forming  an  image  based  on  this  two-dimensional  matrix.  Formation  of  a  correlation 
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image  is  more  complicated,  but  the  basic  approach  is  the  same:  establish  a  time  gate,  calculate 
the  correlation  coefficient  between  adjacent  A-scan  in  the  three-dimensional  matrix,  write  these 
correlations  into  a  two-dimensional  matrix,  and  form  a  correlation  image  based  on  these 
correlations.  As  will  be  discussed,  complications  arise  from:  initial  alignment  of  A-scans,  local 
alignment  of  signals,  skipping  signal  measurement  positions,  and  image  formation  using 
matrices  with  blank  elements. 

The  correlation  approach  addressed  here  contains  at  least  an  element  of  matched  filtering. 
Matched  filtering  is  a  technique  commonly  used  in  radar  to  optimize  the  signal  to  noise  ratio 
(SNR).  Matched  filter  theory  states  that  when  choosing  a  finite  impulse  response  filter  to 
optimize  SNR  of  a  received  signal,  the  ideal  filter  will  have  the  same  frequency  content  as  the 
defect  signal  to  be  detected.  The  application  of  matched  filtering  in  radar,  optics  and  acoustics 
has  an  extensive  history,  whether  applied  classically  (see,  e.g.,  references  4-11)  or  adaptively 
(see,  e.g.,  references  12-17).  However,  challenges  still  remain  relative  to  a  priori  determination 
of  the  frequency  content  of  unknown,  and  yet  to  be  detected,  defect  signals.  As  applied  to  defect 
detection  using  ultrasonics,  implementation  of  the  matched  filter  is  comparable  to  cross- 
correlating  the  expected  defect  signal  (the  “template”)  with  each  received  A-scans  and  using 
high  correlation  values  as  a  basis  for  detection.  The  efficacy  of  this  technique  as  applied  to 
ultrasonics  is  directly  related  to  the  quality  of  the  template,  that  is,  the  filter’s  efficacy  is  related 
to  the  operator’s  ability  to  establish  of  the  frequency  content  of  the  defect  signal  to  be  detected. 
Accounting  for  the  influence  of  scattering  by  an  unknown  defect  on  the  frequency  content  of  the 
defect  signal  is  particularly  difficult.  In  the  correlation  approach  presented  in  this  paper,  each 
signal  is  correlated  with  its  nearest  neighbors.  In  essence,  each  gated  A-scan  is  taken  to  be  a 
template  for  application  of  matched  filtering  to  gated  A-scans  measured  at  neighboring 
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measurement  positions,  with  high  correlations  again  used  as  a  basis  for  detection.  To  some 
degree,  the  neural  processes  of  a  human  observer  are  being  mimicked  by  the  correlation 
approach.  For  example,  an  observer  of  a  B-scan  image  will  notice  local  similarities  in  A-scans 
making  up  the  B-scan  -  even  if  distorted  in  unpredictable  or  unknown  ways,  the  observer  will 
detect  the  similarity. 

The  remainder  of  this  paper  deals  with  the  details  of  the  correlation  approach.  The  paper 
begins  by  addressing  the  procedures  used  in  measuring  A-scans  and  for  creating  simulated  A- 
scans.  Correlation  calculations  and  issues  associated  with  signal  alignment  and  correlation 
image  display  are  then  addressed.  Correlation  images  are  presented  in  the  results  section  based 
on  the  measured  and  simulated  A-scans.  C-scan  images  are  also  presented  for  comparison 
purposes  with  the  new  correlation  images.  The  paper  closes  with  a  summary  and  discussion 
section. 

II.  Measured  and  simulated  A-scans 

A-scans  were  measured  in  an  immersion  mode  from  three  samples:  a  stainless  steel 
plate;  a  stainless  steel  annular  ring  with  nominal  inside  and  outside  diameters  of  40  mm  and  70 
mm,  respectively;  and  an  aluminum  sample  fabricated  to  represent  an  aircraft  structure.  A 
Panametrics  10  MHz,  14”  diameter,  4”  focal  length  transducer  was  used  in  scanning  the  plate 
sample  at  normal  incidence  on  a  45  x  62  regular  grid  with  0.5  mm  between  measurement 
positions.  Data  from  the  plate  sample  was  used  to  show  grain  noise  correlation  images  at 
different  measurement  position  spacing  and  as  a  basis  for  creating  simulated  noise  signals.  The 
annular  ring  sample  was  machined  to  add  five  equal-depth  flat  bottom  holes.  A  Panametrics  1 5 
MHz,  14”  diameter,  6”  focal  length  transducer  was  used  for  a  polar  scan  at  normal  incidence  on 
the  annual  sample  with  2  degrees  between  radial  scan  lines  and  0.762  mm  between 
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circumferential  scan  lines.  Data  from  the  annular  sample  was  used  to  show  polar  correlation 
images  and  to  present  correlation  images  for  flat-bottom-hole  signals  within  grain  noise.  The 
aluminum  sample  was  comprised  of  two  aluminum  plates  fastened  together  with  EDM  notches 
of  various  lengths  extending  from  selected  fastener  holes.  For  the  aluminum  sample,  oblique 
incident  pitch-catch  measurements  were  made  using  a  pair  of  Panametrics  10  MHz,  V2”  diameter, 
transducers,  one  flat  and  one  with  a  focal  length  of  3”.  Transducers  were  arranged  5  degrees 
from  normal.  The  aluminum  sample  with  simulated  cracks  was  used  to  demonstrate  the  use  of 
correlation  images  for  the  detection  of  cracks  based  on  lack  of  correlation  in  transmitted  signals. 

A-scans  were  also  created  using  simulated  defect  signals  buried  within  simulated  grain 
noise.  Gaussian  damped  sinusoids  were  used  to  simulate  defect  signals  with  the  appropriate 
frequency  content.  To  simulate  planar  defects  with  varying  aspect  ratio  and  at  various  signal-to- 
noise  ratios,  defect  signals  were  scaled  over  certain  spatial  regions  and  added  to  the  simulated 
noise  signals.  This  process  is  represented  in  equation  form  for  a  signal  measured  at  the  ij 
position  in  a  simulated  raster  scan  as  follows: 

x(i,  j,  t)  =  n(i,  j,t)+  B(i,  j)y(t)  ( 1 ) 

where  n(i,  j,  t)  represents  simulated  acoustic  noise  before  the  addition  of  a  defect  signal,  y{t) , 
scaled  by  B(i,j),  a  scale  factor  which  can  be  varied  over  lateral  space,  and  x(i,j,t )  is  the  result, 
representing  the  simulated  measurement  data.  The  computation  is  carried  out  over  all  noise 
signals.  Simulated  grain  noise  signals  were  generated  to  mimic  the  maximum  extreme  value 
(gated  peak  value)  distribution  and  correlation  coefficient  distribution  associated  with  the 
measured  noise  from  the  stainless  steel  plate  sample.  [18]  The  basic  procedure  involves  creating 
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each  simulated  A-scan  as  a  scaled  version  of  previously  generated  signals  plus  a  new  random 
noise  component.  Desired  correlation  values  between  signals  are  drawn  at  random  from 
correlation  coefficient  distributions  established  using  the  measured  noise.  The  maximum 
extreme  value  distribution  associated  with  the  simulated  signals  is  controlled  through  the 
standard  deviation  of  the  random  number  generator,  by  scaling  the  signals  to  control  signal 
energy,  and  by  filtering  the  signals  to  yield  the  same  frequency  content  as  the  measured  signals. 
The  details  of  this  procedure  are  beyond  the  scope  of  this  paper,  but  will  be  documented 
elsewhere. 

III.  Formation  of  correlation  images 
A.  Basic  correlation  coefficient  calculation 

The  correlation  of  interest  is  the  spatial  cross-correlation  calculated  between  gated  A- 
scans  measured  at  adjacent  measurement  positions.  Consider  an  Ax  M  scan  with  each  A-scan  T 
points  long  and  using  the  matrix  x  =  x(i,j,t )  i  =  1...A  j  =  \...M  t  =  0...T  - 1  to  hold  the  A- 
scans  (see  Fig.  1).  Calculation  of  the  correlation  coefficient  between  A-scans  is  then  given  by 
the  following  equation  where  p  =  p{i,j,8r,8c,r )  is  a  sample  estimate  of  the  associated  expected 
value:  [19] 


X  Mb  j >  {-T)-mx  (*>  ./)IX(Z  +  8c,j  +  8r,t)-mx  (i  +  Sc,j  +  Sr )] 


t=tj 

,  ZWc/T- 

V=h 

-r)-mx{Uj)f  ^ 

1  (f 

I  Xtx(z  +  +  Sr,t)~  mx(i  +  Sc,j  +  8r  )]2 

t=tt 

(2) 


6 


Normalization  removes  the  scale  dependence  and  restricts  p  to  the  range  -1  to  1.  In  Eq.  (2),  the 
summation  range  in  the  time-domain  defines  the  portion  of  the  signal  (the  time  window  or  gate) 
of  interest,  mx  and  my  are  mean  values  calculated  over  the  gate,  r  controls  the  lag  or  temporal 

shift  between  the  two  signals,  and  8  is  a  spatial  shift  parameter.  Throughout  the  paper,  t  is  used 
as  a  discrete  index  referring  to  the  temporal  direction.  With  8r  =1,  8 c  =  0 ,  row  correlation  are 

calculated  between  adjacent  A-scan,  that  is,  A-scans  measured  at  the  jth  and  jth  + 1  positions  in 
the  ith  scan  row.  Similarly,  with  8r  =  0,  8C  =  1 ,  column  correlations  can  be  calculated. 

Correlations  for  all  possible  adjacent  signal  combinations  in  an  N  x  M  raster  scan  can  be 
established  using  a  computation  loop  over  i  and  j  with  the  spatial  shift  applied  sequentially  to  i 
and  j.  Correlations  can  also  be  calculated  with  signals  systematically  skipped  by  setting  the 
spatial  shift  parameters  to  an  integer  greater  than  one,  for  example,  with  8r  =2  8 c  =  0 ,  row 

correlations  are  calculated  for  every  other  signal. 

B.  Signal  alignment 

As  a  point  of  motivation,  consider  the  same  material  solid-state  welding  problem  where 
the  front  surface  of  the  part  is  nominally  flat  and  parallel  to  the  weld  plane. [20,21]  Signal 
misalignment  can  result  from  errors  in  the  mechanical  scan,  systematic  or  random  variations  in 
the  front  surface,  and  systematic  or  random  weld-surface  irregularities.  As  such,  signal 
alignment  will  be  addressed  in  terms  of  global  alignment  using  front  surface  reflections  and  then 
local  alignment  based  on  gated  A-scans. 

Global  signal  alignment  for  backscatter  is  done  by  aligning  the  leading  edge  of  each  front 
surface  reflection.  In  general,  when  interrogating  the  volume  of  a  component  for  defects,  the 
front  surface  reflection  is  blown  off  scale,  making  most  of  the  front  surface  reflection  truncated 
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and  difficult  to  use  in  a  signal  processing  sense.  Nevertheless,  the  signals  can  be  aligned  based 
on  the  front  surface  reflection  leading  edge  by  applying  any  number  of  signal  processing 
approaches,  including  alignment  of  trigger  points  as  used  here.  Global  alignment  yields  B-scan 
images  aligned  at  the  front  surface  reflection,  providing  a  basis  for  evaluating  the  remaining 
systematic  or  random  misalignment  in  the  weld  plane  signal.  This  global  alignment  is  only 
achieved  easily  and  unambiguously  under  certain  measurement  configurations  and  sample 
geometries. 

Local  signal  alignment  can  be  used  to  form  locally  aligned  B-scan  images  or  to  form 
correlation  coefficient  images  based  on  maximum  correlations.  In  either  case,  the  core  process 
involves  forming  a  cross-correlelogram  between  two  signals.  The  cross-correlelogram,  which 
looks  at  the  correlation  between  signals  as  one  signal  is  shifted  in  time  relative  to  the  other 
signal,  is  given  by  a  plot  of  p(r)  versus  r .  From  this  plot,  the  maximum  correlation  and  the 
shift  to  the  maximum  correlation  can  be  extracted. 

To  perform  local  signal  alignment  for  B-scan  images,  each  gated  signal  is  aligned  to  a 
single  signal  which  has  been  chosen  as  the  reference  alignment  signal.  Any  signal  in  a  scan 
could  be  used  as  the  reference  signal;  however,  we  typically  choose  the  signal  that  has  the 
greatest  voltage  within  the  time  gate.  Alignment  is  achieved  based  on  the  time-shift  to  the 
maximum  correlation  between  each  signal  and  the  reference  signal  as  discussed  above.  Again, 
consider  same  material  solid  state  welding  between  two  polycrystalline  metals.  If  the  weld  is 
inadequate  resulting  in  an  echo  signal  from  the  weld  plane,  the  alignment  procedure  will  align 
the  weld  plan  signals  and  yield  B-scan  images  with  a  relatively  smooth  stripe  at  the  time 
corresponding  to  the  weld  plane  position.  Conversely,  a  perfect  weld  will  yield  only 
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backscattered  grain  noise,  and  the  alignment  procedure  will  result  in  B-scans  which  show  at  most 
short  stripes  randomly  distributed  in  time. 

To  perform  local  signal  alignment  to  be  used  in  correlation  images,  we  calculate  the 
cross-correlelogram  between  each  gated  signal  and  its  nearest  row  and  column  neighbors.  The 
maximum  correlation  value  is  then  extracted  from  the  plot  of  p(r)  versus  r  for  each  pair  of 
signals.  For  an  N  x  M  scan  where  correlations  are  calculated  with  no  signals  skipped,  there  will 
be  N(M  - 1)  row  correlations  and  (N  - 1  )M  column  correlations,  for  a  total  number  of 
2NM-M-N  correlations  which  will  be  used  as  a  basis  for  forming  the  correlation  image. 
Figure  2  shows  a  schematic  representation  of  a  4  x  4  scan  (left  figure)  with  open  circles  at  the  16 
measurement  positions.  The  R’s  and  C’s  are  used  to  represent  the  positions  conceptually 
assigned  to  the  correlation  between  adjacent  row  signals  and  adjacent  column  signals, 
respectively.  Also  shown  in  the  figure  is  a  portion  of  a  polar  scan  (right  figure)  with  the  R’s  and 
C’s  now  representing  the  positions  assigned  to  radial  and  circumferential  correlations, 
respectively,  and  the  open  circles  again  at  the  measurement  positions. 

C.  Image  display  issues 

Whether  before  or  after  signal  alignment,  there  are  issues  associated  with  displaying  the 
N{M  - 1)  +{N  -\)M  correlation  values  as  an  image.  For  discussion  purposes,  again  consider  a 
4x4  raster  scan.  Nearest  neighbor  correlations  will  result  in  a  4  x  3  set  of  row  correlations  and  a 
3  x  4  set  of  column  correlations  as  is  apparent  from  Fig.  2.  At  this  point,  the  easiest  approach 
would  be  to  make  two  correlation  images:  one  for  the  row  correlations  and  one  for  the  column 
correlations.  Similarly,  radial  and  circumferential  correlation  images  could  be  made  for  a  polar 
scan.  This  procedure  yields  two  correlation  images  which  allow  separate  evaluation  of  row  and 
column  correlations. 
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The  issues  associated  with  forming  a  single  correlation  image  can  be  addressed  by  first 
thinking  in  terms  of  writing  the  row  and  column  correlations  into  a  single  matrix.  For  the  4x4 
example,  to  maintain  spatial  registration,  the  4  x  3  row  correlations  and  the  3  x  4  column 
correlations  would  be  written  into  a  7  x  7  matrix  (see  Fig.  2).  For  this  example,  the  difficulty  is 
that  only  24  correlation  coefficients  are  written  into  a  49  position  matrix,  leaving  25  blank 
elements.  These  blank  elements  spatially  correspond  to  the  measurement  positions  (open  circles 
in  Fig.  2)  and  to  the  centroid  of  the  region  defined  by  each  set  of  four  measurement  positions 
(solid  circles  in  Fig.  2).  In  general  terms,  the  N(M  - 1)  +(n-\)m  row  and  column  correlations 
are  written  into  a  (2 N  —  l)  x  (2 M  - 1)  matrix,  leaving  2 NM  —  N  —  M  + 1  blank  elements.  An 
image  based  directly  on  a  matrix  with  blank  elements,  with  each  blank  element  assigned  the 
same  value,  is  visually  unsatisfying  with  useful  correlation  values  displayed  within  a  regular 
array  of  mono-colored  pixels. 

A  variety  of  solutions  exist  to  this  image  display  problem.  One  potential  solution  would 
be  to  assign  each  blank  element  an  average  of  nearest  neighbor  values.  The  resultant  image 
would  be  composed  of  pixels  associated  with  row  correlations  and  column  correlations  and 
pixels  which  are  an  average  of  row  and  column  correlations.  A  second  approach  is  to  shift  the 
column  correlations  up  one  row.  This  approach  compromises  spatial  registration;  however,  as 
shown  in  the  results  section,  the  resultant  images  are  useful  for  qualitative  evaluation.  A  third 
approach  which  preserves  spatial  registration  is  to  consider  each  correlation  value  as  being 
centered  in  a  diamond-shape  region.  This  approach  can  be  applied  to  both  Cartesian  and  polar 
scan  data  as  depicted  in  Fig.  2.  For  the  4  x  4  example,  the  image  would  have  4  rows  of  3 
diamonds  (row  correlations)  interleaved  with  3  rows  of  4  diamonds  (column  correlations).  To 
facilitate  display  of  these  diamond  regions,  a  matrix  with  a  much  higher  number  of  elements  than 
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the  original  correlation  matrix  is  used  to  approximate  the  diamond-shaped  regions,  that  is,  the 
diamonds  are  discretized  (see  Fig.  3).  Each  element  in  the  high  density  matrix  is  assigned  a 
value  based  on  which  diamond  that  element  falls  within.  The  pixels  outside  of  the  diamond- 
filled  regions  are  set  to  a  single  (arbitrary)  value.  As  shown  in  the  blown-up  polar  correlation 
images  in  Fig.  4,  a  close  look  at  each  diamond  region  in  a  correlation  image  reveals  this 
discretization  with  each  diagonal  line  having  a  saw-toothed  appearance  at  the  scale  of  the  high 
density  matrix  and  the  outside  edges  of  the  image  having  a  saw-toothed  appearance  at  the  scale 
of  the  diamond  regions. 

IV.  Results 

Example  results  were  chosen  to  cover  measured  and  simulated  data,  xy  and  polar  scans, 
planar  and  localized  defects,  and  detection  based  on  increased  and  decreased  correlation, 
respectively.  In  each  case,  companion  C-scan  images  are  shown  for  comparison  purposes. 

Figure  5  shows  a  correlation  image  (top)  and  C-scan  image  (bottom)  for  the  stainless 
steel  annular  ring.  Data  were  measured  in  a  polar  scan  with  the  correlation  image  displayed 
using  the  diamond  region  format.  The  regions  of  higher  correlation  and  higher  amplitude  near 
the  inside  and  outside  diameter  are  due  to  multiple  reflected  signals  from  surfaces  not  in  the 
plane  of  the  flat-bottom-holes.  The  flat-bottom-hole  diameters  starting  with  the  top  hole  and 
moving  clockwise  are  1/16”,  2/16”,  3/16”,  3/16”,  and  3/16”.  In  both  images,  the  five  flat- 
bottom-holes  are  clearly  visible.  Detection  in  the  correlation  images  is  possible  due  to  increased 
correlation  between  signals  containing  echo  signals  from  a  given  flat-bottom-hole. 

Figure  6  shows  three  images  based  on  simulated  data  generated  as  described  above.  The 
top  figure  is  an  image  of  the  scale  factor  /?(/,/) .  This  image  represents  a  SNR  mask  which 
reveals  the  location  and  SNR  of  each  of  the  simulated  planar  defects  which  were  superimposed 
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on  simulated  acoustic  noise.  SNR  is  defined  here  as  the  ratio  of  the  maximum  extreme  value  of 


the  signal  to  be  added  to  the  noise  and  the  expected  maximum  extreme  value  of  the  noise.  The 
SNR  for  each  defect  is  indicated  in  the  figure  caption.  Notice  that  the  highest  SNR  is  less  than 
one.  The  middle  image  in  Fig.  6  is  a  C-scan  image.  Even  at  the  very  low  SNR’s,  given  the 
lateral  extent  of  the  defects,  three  out  of  five  defects  are  visible  in  the  C-scan  image.  A  fourth 
defect  also  visible,  especially  given  the  bias  created  by  knowing  the  defect  locations.  The  lower 
image  is  a  correlation  image  formed  using  the  interleaving  approach.  In  this  image,  four  of  the 
defects  are  easily  detected  with  the  region  around  the  fifth  defect  also  arguably  of  concern. 

The  final  example  is  based  on  measured  oblique  incidence  pitch-catch  data  from  the 
aluminum  sample.  The  schematic  in  Fig.  7  shows  the  sample  geometry.  The  sample  is  made  up 
of  two  aluminum  plates  joined  together  with  two  rows  of  eight  fasteners  with  sealant  between  the 
two  plates.  EDM  slots  were  machined  to  extend  from  four  of  the  fastener  holes  in  the  outer 
aluminum  plate  (the  skin)  at  four  different  lengths  as  shown  in  the  figure.  Images  are  based  on 
the  back  surface  reflection  (signal  C)  from  the  lower  plate  (the  substructure).  Similar  results  are 
obtained  by  monitoring  signal  B  from  the  interface  (the  faying  layer).  In  this  case,  interaction  of 
the  acoustic  beam  with  the  EDM  slots  causes  amplitude  reduction  and  distortion  of  the 
transmitted  signal.  The  reduction  in  amplitude  makes  the  presence  of  the  cracks  apparent  in  the 
C-scan  image.  The  distortion  leads  to  a  reduction  in  signal  correlation,  making  the  cracks 
detectable  in  the  correlation  image.  The  three  longer  slots  are  easily  detected  and  sized  in  either 
image  with  the  1  mm  crack  difficult  to  detect  without  optimized  measurements  and  additional 
processing. 
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V.  Summary  and  discussion 

The  correlation  approach  for  defect  detection  outlined  in  this  paper  is  attractive  as  a 
potential  compliment  to  the  classical  C-scan  approach  for  a  number  of  reasons  including:  the 
correlation  approach  is  scale  independent;  the  correlation  approach  relies  on  signal  shape,  while 
the  C-scan  approach  uses  only  the  peak  value;  correlation  images  are  calculated  based  on  the 
same  A-scans  used  in  forming  C-scan  images;  correlation  images  have  a  qualitatively  similar 
appearance  to  C-scan  images,  which  facilitates  inspector  acceptance;  and  at  very  low  SNR,  the 
correlation  approach  shows  the  potential  to  outperform  the  C-scan  approach  for  detecting  defects 
that  show  some  lateral  extent. 

In  terms  of  additional  development  on  the  correlation  approach,  work  is  underway  on  a 
number  of  fronts.  A  formal  ROC  analysis,  which  is  underway,  should  shed  light  on  the 
performance  of  the  correlation  approach  alone  and  in  combination  with  the  C-scan  approach 
relative  to  the  performance  of  the  C-scan  approach  alone.  A  matched  fdter  approach  will  also  be 
included  in  the  analysis.  The  correlation  approach  is  amplitude  independent  but  is  not 
measurement  system  and  material  independent.  Preliminary  assessment  has  shown  very  weak 
dependence  of  correlation  values  on  the  beam  field  near  the  focal  region  of  a  circularly  focused 
piston-source  transducer.  Future  work  will  include  a  more  in  depth  study  of  the  influence  of 
beam  field  variations,  frequency  content,  and  material  morphology  on  the  correlation  approach. 
Initial  work  has  also  been  done  in  the  area  of  statistical  detection  applied  to  the  correlation 
approach  with  the  goal  of  creating  threshold  images  based  on  a  specified  false-call  rate.  We 
conclude  by  noting  that  this  is  only  a  partial  list  of  issues  to  be  considered  in  developing  and 
evaluating  this  correlation-based  approach  to  image  formation  and  defect  detection. 
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Figure  1.  Schematic  representation  of  a  typical  raster  scan  and  data  matrix. 
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Figure  2.  Schematic  representation  of  Cartesian  and  polar  scans.  Open  circles  represent 
measurement  positions.  R’s  and  C’s  represent  the  conceptual  positions  of  row  and  column 
correlations  for  the  Cartesian  scan  and  radial  and  circumferential  correlations  for  the  polar  scan. 
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Figure  3.  Correlation  image  display  approach  showing  the  highly  discretized  approach  used  to 
represent  diamond  regions  surrounding  row  and  column  correlations. 
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Figure  4.  Example  correlation  images  using  the  diamond  region  approach  for  a  polar  scan. 
Each  diamond  region  represents  either  a  radial  or  circumferential  correlation  coefficient.  Images 
are  shown  at  high  magnification  in  order  to  reveal  pixelization  details. 
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Figure  5.  Example  correlation  image  (top)  and  C-scan  image  (bottom)  showing  5  flat  bottom 
holes  in  a  stainless  steel  annular  ring.  The  correlation  image  was  formed  using  the  diamond 
region  approach.  The  darker  regions  in  the  images  represent  high  correlations  or  high 
amplitudes. 
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Figure  6.  C-scan  image  (middle)  and  correlation  image  (bottom)  showing  detection  of  low  SNR 
simulated  defects  in  computer  generated  acoustic  noise.  The  darker  regions  in  the  images 
represent  high  correlations  or  high  amplitudes.  The  position  and  SNR  of  each  defect  is  shown  in 
the  upper  image.  The  correlation  image  was  formed  using  the  interleaving  approach. 
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Figure  7.  Example  of  crack  detection  in  pitch-catch  based  on  decreased  correlations.  Darker 
regions  in  the  C-scan  image  represent  decreased  amplitude.  Lighter  regions  in  the  correlation 
image  represent  decreased  correlations.  The  correlation  image  was  formed  using  the  interleaving 
approach. 
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